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Abstract 

The Large Area Telescope (LAT) on-board the Fermi Gamma-ray Space Tele- 
scope began its on-orbit operations on June 23, 2008. Calibrations, defined in a generic 
sense, correspond to synchronization of trigger signals, optimization of delays for latch- 
ing data, determination of detector thresholds, gains and responses, evaluation of the 
perimeter of the South Atlantic Anomaly (SAA), measurements of live time, of abso- 
lute time, and internal and spacecraft boresight alignments. Here we describe on-orbit 
calibration results obtained using known astrophysical sources, galactic cosmic rays, 
and charge injection into the front-end electronics of each detector. Instrument re- 
sponse functions will be described in a separate publication. This paper demonstrates 
the stability of calibrations and describes minor changes observed since launch. These 
results have been used to calibrate the LAT datasets to be publicly released in August 
2009. 
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1 Introduction 

The Fermi Gamma-ray Space Telescope, hereafter Fermi, represents the next generation of 
satellite-based high-energy gamma-ray observatory. The Fermi satellite hosts two instru- 
ments: the Large Area Telescope (LAT) pj and the Gamma-ray Burst Monitor (GBM) [2]. 
The former employs a pair-conversion technique to measure photons from 20 MeV to energies 
greater than 300 GeV, while the latter uses Nal and EGO scintillation counters to record 
transient phenomena in the sky in the energy range from 8 keV to 40 MeV. The LAT has no 
consumables, and a very stable response unlike its predecessor, the Energetic Gamma Ray 
Emission Telescope (EGRET) [3j. For the energy range above 10 GeV the sensitivity of the 
LAT is at least one order of magnitude greater than that of EGRET, allowing the sky to be 
explored at these energies essentially for the first time [1]. 

The LAT consists of a tracker /converter (TKR) for direction measurements [HEIEIE], 
followed by a calorimeter (GAL) for energy measurements [S]. Sixteen TKR and GAL 
modules are combined to form sixteen towers, which are assembled in a 4x4 mechanical 
support structure. An anticoincidence detector (AGD), enclosed by a micrometeoroid shield, 
surrounds the TKRs and rejects charged cosmic-ray background [9l[T0]. The LAT has about 
one million detector readout channels. 

On-orbit calibrations relate to all aspects of LAT measurements and data analysis results, 
from absolute timing to energy and direction measurements for individual events, to fluxes 
and positions of gamma-ray sources. 

The accurate timestamps of the LAT are obtained using the Global Positioning System 
(GPS) of the Fermi spacecraft, which provides timing and position information. Those are 
needed for phase folding pulsars and correlating gamma-ray observations with those at other 
wavelengths. 

As the number of photons in an observation increases, the centroid of their spatial dis- 
tribution becomes better measured, and eventually the error is dominated by uncertainties 
in the alignment of the LAT, both internal and with respect to the Fermi spacecraft. 

Source localization at GeV energies enables the LAT to resolve bright, adjacent sources 
previously labeled as unidentified [3] and will help elucidate the origin of gamma-ray emis- 
sions from galactic cosmic rays accelerated in supernova remnants. 

The energy calibrations at higher energies are of utmost importance for detection of dark 
matter particle signals. Some extensions to the Standard Model predict narrow spectral 
lines due to the annihilations of as-yet unknown massive particles. For detecting and char- 
acterizing these features, accurate energy determination is vital. Even though the LAT was 
designed to measure gamma-rays it can also study, though not separately, the cosmic-ray 
electron and positron spectra. Gosmic-ray electron and positron spectra and intensities may 
also contain signatures for new physics. In this case, energy calibrations and position deter- 
mination using extrapolated tracks into the GAL play an important role. At lower energies, 
broad features in the photon energy spectra of active galactic nuclei, or supernovae remnants 
originating from pion decays and bremsstrahlung may help unravel outstanding questions 



concerning particle acceleration in these sources. 

The purpose of this paper is to document the on-orbit calibration procedures used by the 
LAT; it begins with an overview of calibrations in Section [2j Details on trigger, ACD, CAL 
and TKR cahbrations are described in Sections[3llll[5]and[6l respectively. The evaluation and 
updates to the perimeter of the South Atlantic Anomaly (SAA) are presented in Section [7] 
and measurements of live time are discussed in Section [Si The results from absolute timing 
follow in Section O Finally, internal and spacecraft boresight alignments are explained in 
Section [TOl We conclude with a table that summarizes the calibration results in Section [TTl 
Assessment of the current LAT performance is described in a separate publication [11], while 
tests performed at particle accelerators are presented elsewhere [12], [131 El- The section 
describing each calibration in shown in the last column of Table [H 

2 Overview of LAT calibrations 

For this paper, the word calibration represents synchronization of trigger signals, optimiza- 
tion of delays for latching data, determination of detector thresholds, gains and responses, 
evaluation of the perimeter of the South Atlantic Anomaly (SAA), measurements of live 
time and of absolute time and internal and spacecraft boresight alignments. 

Table [1] summarizes all calibration types, classified by category. Configuration refers to 
operational settings that define the state of the hardware, which are fixed before data are 
acquired. Calibrations are related to quantities that can change after data are processed and 
analyzed. As shown in Table [U for some types of calibrations, updates are not as frequent 
as the acquisition of the relevant data. Operationally, calibration data are acquired in two 
distinct modes (see Section [3TT1) : 

1. Dedicated, meaning that the trigger, detector and software filter settings are incom- 
patible with nominal science data taking. 

2. Continuous, meaning the trigger and software can, with only a small penalty in live 
time, acquire specialized data that is used to calibrate, or, more generally, monitor the 
performance of the LAT during nominal science data-taking. 

Within a run the LAT acquires data with fixed instrument configurations. For the most part, 
the dedicated calibration runs are concerned with characterizing the electronics' response to 
known stimuli while the continuous calibrations are aimed at calibrating the electronics 
with a known physics input. The stability of calibrations has been such that operations in 
dedicated-mode amount to approximately 2.5 hours every three months. 

Sea-level cosmic ray muons were used to calibrate the low-energy scales and trigger 
thresholds, but muons do not deposit enough energy in the detector elements to calibrate 
the high-energy scales and high-energy trigger thresholds. Instead, we used charge injection 
into the front-end electronics to calibrate the high-energy scales. Because of rise-time slewing 
effects, the optimal synchronization of trigger signals and optimal delays for data latching are 
energy dependent. We used pre-launch tests to provide a best approximation of the optimal 
trigger timing, and verified and corrected the synchronization and delays with on-orbit data. 
All these calibrations are revisited with on-orbit data selected from galactic cosmic rays and. 
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Table 1: List of LAT calibrations (calib) and configurations (config) that can impact LAT 
scientific results. Note that pedestals are used both as a configuration on-board and as a 
calibration on the ground. Detailed descriptions are given in the sections listed in the last 
column. Frequency of updates correspond to current best estimates for the period beyond 
the first year of operations. 



as it will be shown later, there are only minor deviations when comparing results prior and 
after launch. 

3 Overview of trigger and readout 

The trigger and data readout system controls the composition and flow of data from the 
source in the detector elements to the Solid State Recorder (SSR) of the Fermi spacecraft. 
There are two distinct stages: the latching and movement of data from the front-end detector 
elements to a LAT Event Processing Unit (EPU), and the processing of the data on the EPU 
and its subsequent transfer to the SSR. The first stage is controlled by hardware, namely the 
trigger, and the second stage by software. A full description of the LAT multilevel trigger 
and data readout can be found elsewhere [H 

3.1 Acquisition modes: dedicated or continuous 

During dedicated calibration runs the trigger system and detector electronics are configured 
to acquire data useful for calibrating the thresholds and responses of detectors and synchro- 
nizing the arrival times of signals from various parts of the detectors. We call these runs 
dedicated because the trigger and detector configurations necessary to acquire these special- 
ized data are incompatible with acquiring our primary science data (i.e. high-energy gamma 
rays from celestial sources). During continuous calibration acquisitions, the trigger system 
and detector electronics are configured to detect and latch not only those events thought to 
be gamma rays, but also those useful for calibration. 

As used on orbit, the LAT trigger system takes inputs from the ACD, TKR, and CAL 
front-end electronics and from a programmable, internal periodic trigger. Under control of 
LAT flight software, the periodic trigger can be used with a programmable charge- injection 
system to calibrate the detector electronics, or it can be used simply to read out the detector 
front-ends at a specified cadence. 

The programmed control under the charge injection system is used only in dedicated 
calibration runs. There the trigger is configured to collect a specified number of events 
at a particular rate, regardless of input from the detectors. In most cases, the trigger is 
configured to instruct the detector electronics to inject a known, programmable amount of 
charge at a specified time relative to each trigger. In other cases, no charge is injected. 
For each sequence of these events, flight software sets configuration registers throughout the 
instrument that control, for example, the amount of charge injected. By collecting a number 
of events with the same injected charge, and then varying the amount of injected charge, 
the electronic response can be accurately calibrated. 

Other dedicated calibration runs and all continuous calibration runs have the instrument 
triggered by signals from the detectors, not a programmed sequence, and charge injection is 
not used. For example, dedicated calibration runs in this mode are those used to synchronize 
trigger signals from the detectors. Data are collected using a sequence of configurations that 
sweep values through the various registers that control the delay times of the trigger signals. 
The data captured are then analyzed on the ground to identify the delay settings that best 
synchronize the arrival times of the trigger signals from each detector system. 



During continuous calibration runs, the trigger is configured to detect and latch both 
events thought to be gamma rays and other events useful for calibration. Based on the 
signals recieved from the detectors, the trigger specifies the type of ACD and CAL readout. 
There are four types (see Table [2]) formed by enabling or disabling zero-suppression, in which 
values are committed to the data stream only if they exceed a programmable threshold, and 
selecting CAL single-range or four-range mode, in which the electronics takes respectively 
only the "best" range or all four values (see Section EH for details). The ACD does not 
support a readout mode to select both readout ranges. It does, however, respect the zero- 
suppress mode. 



Purpose 


Readout Mode 


Select photons 


Zero-suppressed 
and best range 


Cross-calibration of the energy ranges (CAL) 


Zero-suppressed 
and four-range 


Calibration of the MIP peak (CAL,ACD) 


Zero-suppressed 
and best range 


Monitor pedestals (CAL, ACD) and noise occupancy (TKR) 


Non-zero-suppressed 
and four-range 



Table 2: Readout modes used for continuous calibrations (nominal science operations). 



After an event passes the hardware trigger it is inspected by on-board software filters, 
each configured to identify events likely to be useful for one or more scientific or calibration 
purposes. If any filter accepts an event, it is included in the LAT data stream and forwarded 
to the SSR for transmission to the ground. In addition, the filters can be configured to 
allow a fraction of events that would otherwise have been rejected to be included in the data 
stream. Such events are used primarily to study filter performance. The input trigger rate 
of about 2.2 kHz, averaged over many orbits, is reduced to about 450 Hz by the software 
filters. Table [3] describes filter types, purpose and average output rate. 

3.2 Synchronization of trigger signals and delays for latching data 

Fast trigger signals (~ few hundred ns) from the detectors must be synchronized with respect 
to each other for the trigger to operate efficiently. There are five such signals: veto and high 
level discriminator in the ACD (see SectionH]), low and high-energy in the CAL (see Section[5]) 
and TKR (see Sec|6]). The earliest arriving trigger signal initiates a readout. The arrival 
times of the other trigger signals relative to this signal are captured in the event data with 50 
ns precision (i.e. the period of the LAT system clock), allowing a direct comparison of trigger 
signals on an event-by-event basis. Once the relative timing is determined, the settings of 
various delay registers (i.e. the instrument configuration) are modified to synchronize these 
trigger signals in subsequent data acquisitions. There are two types of delay associated with 
the trigger: the fast trigger signal delay and the delay for latching the data. To synchronize 



Filter Type 


Purpose 


Average Rate (Hz) 


Lramma 


select gamma-ray candidates 


410 


and events > 20 GeV 




Heavy ion 


calibration of high-energy scales, 


2.5 


MIP 


select non-interacting protons 


(nominal) 




10 (dedicated-mode) 


Diagnostic 


filter performance, background 


22 



Table 3: On-board filters used to select events for calibration acquisitions. Rates depend on 
geomagnetic and other orbital variations. Here we list average rates. 



Measurement 


ACD Delay 


TKR Delay 


GAL Delay 


Ground 


800 ns (16 ticks) 


250 ns (5 ticks) 





Orbit 


750 ns (15 ticks) 


200 ns (4 ticks) 






Table 4: Trigger signal delays for measurements on the ground and on orbit. 

the detectors we use the time of arrival of trigger signals, referred to as condition arrival 
times. This is event information containing the number of 50 ns clock ticks that have passed 
between the opening of the coincidence time window and the arrival of a given trigger signal. 
Since the trigger signals for TKR and GAL have significant time walk depending on the 
ratio of signal size to threshold value, it is important to choose an appropriate dataset to 
optimize the timing. We select MIP data for the relative timing of AGD and TKR, since 
the main purpose of the AGD is to reject charged particle background. The relative timing 
between TKR and GAL is optimized using photon candidates selected by the gamma filter. 
The GAL low and high-energy signals are controlled by a single delay and cannot be tuned 
independently. 

Figure [1] shows the arrival times for GAL with respect to the TKR. Negative values 
indicate that the TKR trigger arrived earlier than the GAL trigger signal. Distributions are 
fully contained inside the time coincidence windowlllof 700 ns (14 ticks). The arrival times 
are affected by the ratio of the crystal energy over the threshold, thus influencing the shape 
of the trigger signal curves. The spectrum is slightly altered from its original form due to a 
trigger event selection. There is no difference in arrival times between GAL high- and GAL 
low-energy triggers at energies very far above the two thresholds. 

On the ground, only muons were available for timing calibrations so a small change in 
parameters was observed on orbit. The results from the synchronization of the trigger signals 
are shown in Table |H 

The optimal delay is obtained after analyzing data acquired with fixed delay values. For 
the GAL we fit the MIP peak for each dataset with a fixed delay. Figure [2^ shows an example 
for a delay setting of 750 ns (25 ticks). Figure [2)d shows GAL MIP peak positions for six 



tPrior to launch the time coincidence window was set to 600 ns (12 ticks). 
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Figure 1: On-orbit arrival times of the CAL low-energy (FLE) and high-energy (FHE) trigger 
signals with respect to the TKR trigger. This special dataset was acquired with the time 
coincidence window of 1550 ns (31 ticks). Most of the entries fall within the window used 
for nominal science operations (700 ns or ±13 ticks). 




Figure 2: a) Energy deposited in a single CAL crystal (averaged over all towers) for a fixed 
value of the delay for latching data. The fit (curve) is used to determine the MIP peak 
position in the CAL, b) MIP peak positions for different CAL delay values (all towers). The 
curve is a parabola fit to the data. 



delay values for latching data. The optimal setting is the peak position obtained from the 
fit to the data. There is a data point slightly off the curve due to changes in geomagnetic 
conditions along the orbit. Since data were not recorded at the same location in orbit, the 
MIP selection cuts are designed to keep variations to < 0.1% of the MIP peak value. A 
similar procedure is applied to the ACD. For the TKR, the quantity of interest is not the 
MIP peak position but instead the detector efficiency, which is defined as the number of layer 
hits between the first and the last hit of the track divided by the expected number of layers 
crossed by the track. The optimali^ettings for the on-orbit delays for latching data are 200 
ns (4 ticks), and 2450 ns (49 ticks) for the ACD, CAL and TKR delays, respectively. 



4 ACD calibrations 

The ACD is the LAT first-level defense against the charged particle cosmic-ray background 
that outnumbers the gamma-ray signals by 3-5 orders of magnitude. It consists of 97 separate 
plastic scintillating detectors - 89 scintillator tiles and 8 scintillator ribbons, each viewed by 
two photomultiplier tubes (PMTs) for redundancy. The overall ACD detection efficiency 
is > 0.9997, which is provided by ensuring high uniformity of the detectors' response and 
a large number of photoelectrons. The segmentation is needed in order to minimize pulse 
height variations over the ACD area and to minimize unwanted self- veto due to backsplash of 
soft photons from the developing electromagnetic shower in the CAL. Such self- veto caused a 
significant reduction in effective area in EGRET at high energies [21 [IS]- Detailed information 
about the ACD can be obtained elsewhere P, [TU] . 

ACD calibrations include the determination of the mean values of pedestals, of the sig- 
nal pulse heights produced by single MIP particles in each ACD scintillator, and the veto 
threshold settings; and the high-energy and coherent noise calibrations. All those parameters 
are determined for each PMT, so that each ACD tile or ribbon has two calibrated values. 
Every ACD channel has two ranges, low and high, in order to expand the dynamic range 

^The efficiency for latching TKR data is unchanged up to > 750 ns (15 ticks). 
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Figure 3: a) Pedestal distribution from periodic triggers for a single PMT, b) Long term 
pedestal trending for a single PMT. 

of processed signals. The low range covers the signals below 4-8 MIPs (depending on the 
channel), and the high range extends well above 1000 MIPs. The switching between ranges 
occurs automatically depending on the amplitude of the signal. 

In the future there may be changes to the MIP peak positions due to degradation of 
PMT photocathodes, scintillators or in an optical path between them. A possible way to 
mitigate these effects is to raise the high voltage for the PMTs. Since each high voltage is 
common to groups of 16 or 17 PMTs, this requires re-calibration of all channels belonging 
to that particular PMT group. 

4.1 Pedestals and coherent noise 

Pedestals are offset voltages present at the Analog-to-Digital Converter (ADC) inputs in 
the low and high range readouts. We extract the ACD pedestals for the low range readout 
from the 2 Hz periodic triggers, which currently provide approximately 10,000 random sam- 
ples per orbit. Since a small fraction of these events contains particle signals or electronics 
noise or tails from previous signals in the ACD, we extract the pedestal values by perform- 
ing a Gaussian fit to the central 80% of the pulse height distribution (truncated mean). 
Figure. [HH shows a typical pedestal distribution for a single PMT, where the width of the 
Gaussian (truncated) is about 2-4 pulse height bins (~ 0.01 MIPs). The narrow core of 
the distribution shows the intrinsic electronic noise; tails are residual signals from particles 
near in time to the periodic trigger. The small peak at about 250 is dominated by coherent 
noise contributions. Extracting the pedestal for the high range readout requires a special 
data-taking configuration, which forces a series of randomly triggered events to be read out 
in the high range. Since this configuration is incompatible with regular data taking and the 
pedestals are reasonably stable, these data are only acquired during the quarterly calibration 
periods where the LAT is in dedicated-mode. The data analysis is similar to that of the low 
range pedestals. Figure [3]d shows the long term pedestal trending data for a single PMT. All 
values are plotted relative to the calibration being used, at the time of writing, in the offline 
reconstruction software (Mission week 25). To reduce data volume, we use the low range 
pedestal values on-board. We reject all signals less than 25 counts above pedestal, which 
corresponds to approximately 8 times the electronics noise, or 0.05 times the MIP signal in 
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Figure 4: Readout related pedestal ringing in a single PMT. The truncated mean and RMS 
ADC values obtained using periodic triggers versus the time difference between consecutive 
readout cycles in 50 ns clock ticks. The curve corresponds to a fit of a sinusoidal oscillation 
inside a decaying exponential envelope. 



a typical tile. We also use pedestal values in the flight software data compression algorithm 
(for both ACD and CAL). In this case we reduce the data size by referencing signals against 
pedestal values rather than against zero. 

During ground testing of the ACD we discovered that the readout process causes the 
electronics pedestals to ring. These oscillations occur when elements of its internal circuitry 
resonate at their characteristic frequency. This reproducible effect can be quantified. FigureH] 
shows the difference (A Pedestal) between the coherent noise and the regular pedestal values 
for a single PMT. This truncated mean (see Section l^?Ti) is displayed versus the time difference 
between consecutive readout cycles in 50 ns clock ticks. Data were obtained using periodic 
triggers and the error bars correspond to the RMS ADC values. The curve represents a fit to 
a sinusoidal oscillation inside a decaying exponential envelope, which falls to the electronics 
noise after about 200 ^s. 

As this effect is present in every channel of the ACD, events read out when the pedestal 
peaks at about 50 /is (1000 ticks) after the previous readout can lead to small signals in many 
PMTs. This correction is applied to the offline data to avoid overestimating the total energy 
in the ACD and compromising background rejection and photon selection. For events taken 
at the peak of the coherent oscillation (~ 500 ticks), this calibration reduces the coherent 
noise contribution from 0.05 MIPs per PMT to less that 0.005 MIPs per PMT. The effect is 
so small that we do not need to apply corrections to the on-board processing. Temperature 
dependences in ACD pedestals are negligible. 
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Figure 5: a) Distribution of pedestal subtracted and pathlength corrected signals for a single 
PMT in the ACD; the MIP peak is clearly visible and the curve is a fit to the data, b) Long 
term trending of the MIP peak versus mission week. 

4.2 MIP peak and high range cahbrations 

The MIP peak values are determined using reconstructed tracks pointing at an ACD tile 
or ribbon that recorded a signal. The peak value of the pathlength corrected pulse-height 
distribution for each PMT is the MIP peak for that channel. This calibration is a heuristic 
attempt to quantify the average MIP signal seen in the ACD and not a precise determination 
of all details of the energy deposition in the ACD sensors. Figure [5^ shows a distribution 
of pedestal subtracted and pathlength corrected signals for a single PMT in the ACD. The 
MIP peak at about 600 is clearly seen and the peak below 100 corresponds to soft X- 
ray background. Figure [5]d displays the long term trending of the MIP peak for a single 
PMT, where pedestals from mission week 3 are used as a reference. The stability of the 
measurement is about 10% of the MIP peak for a period of 20 weeks. During the offline 
reconstruction we use the MIP peak calibration values to express raw signal pulse height 
measurements in MIP equivalent values. These are converted into energy using the energy 
deposition in the scintillator (~2 MeV/cm). MIP peaks in the ACD tiles occur between 400 
and 1000 pulse height bins above pedestal, and are determined with an accuracy better than 
5%. 

The high and low range readout are calibrated in an analogous way. The main difference 
is that for the high-range readout we use tracks identified as carbon nuclei by the CAL, since 
the proton MlP-like signals are too small. The algorithm requires the deposited energy in 
the CAL, which is pathlength corrected, to be consistent with that of carbon and the first 
few layers to have similar energy to avoid carbon interacting events. Also, in converting 
from raw pulse heights to MIP equivalent values we allow for non-linearities caused by signal 
saturation in the electronics and scintillators. Figure M shows the deposited energy for the 
full dynamic range of a single PMT. The solid line shows events read out in the low range 
and the dotted line shows events read out in the high range. Figure [7] shows a distribution 
of pedestal subtracted and pathlength corrected signals for carbon in a single PMT. 
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Figure 6: Deposited energy for the full dynamic range of a single PMT. The solid line shows 
events read out in the low range, the dotted line shows events read out in the high range. 
The MIP peak occurs ~ 1.3 MIPs because pathlength corrections are not applied. 




Pathlength corrected signal (ADC) 

Figure 7: Distribution of pedestal subtracted and pathlength corrected carbon signals for a 
single PMT in the ACD in the high range. The carbon peak is clearly visible at ~200 counts 
and the curve shows the fit to the data. 
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Figure 8: Veto turn-on curve for a single PMT. The small number of events (0.1%) < 0.4 
MIPs correspond to overshoots in the ADC value from a previous event. 

4.3 Veto threshold and high-level discriminator 

The on-board thresholds are controlled by DAC settings in the front-end electronics. These 
are calibrated by scanning three settings and measuring the resulting veto threshold in pulse 
height bins. We combine that information with the MIP peak to pulse height scale, or 
carbon peak for the high range, and set the veto and high-level discriminator thresholds as a 
fraction of the MIP and carbon signals, respectively. The veto thresholds for each PMT are 
set separately to an accuracy of about 0.01 MIPs (~ 20 keV) relative to the calibrated MIP 
peak. As shown in FigEl the veto turns on between 0.4-0.5 MIPs and the 50% efficiency point 
is close to 0.45 MIP. For events below < 0.4 MIPs the ADC values can be artificially small, 
because the readout electronics are relatively slow when compared to the veto electronics. 
Events arriving when the electronics are overshooting the return to baseline tend to fire the 
veto discriminator, even though they are not expected to have any ADC counts. The turn-on 
for the high level discriminator occurs between 24-26 MIPs and the RMS width of the carbon 
peak is 20% of the mean value. We have not yet monitored the stability of the carbon peak 
since it needs large statistical samples. At least four months of data are required to obtain 
a carbon peak value with reasonable statistics. 

5 CAL calibrations 

The CAL is designed to measure the energy of incident photons and charged particles, and 
to determine the direction and energy of photons and charged particles for which the TKR 



did not provide direction information, either because their trajectories did not cross the 
TKR or because they did not pair-produce in the TKR. Its imaging properties are also a 
key ingredient in seeding the track reconstruction process in TKR data analysis and in the 
rejection of charged-particle background [1]. 

The CAL consists of 16 identical modules. Each module is composed of 96 CsI(Tl) 
scintillation crystals arranged in a hodoscopic configuration with eight layers each containing 
12 crystals. Each layer is rotated 90° with respect to its neighbors, forming an x-y array. 
Crystals are read out by two dual-PIN-photodiode assemblies, one at each end, that measure 
the scintillation light produced in the crystal. Each photodiode assembly contains a large- 
area photodiode to measure small energy depositions and a small-area photodiode to measure 
large energy depositions. The active areas of the large and small diodes have a ratio of 6:1 
with a spectral response well matched to the scintillation spectrum of CsI(Tl). Each of the 
3072 photodiode assemblies is read out by an amplifier-discriminator ASIC, the GLAST 
Calorimeter Front-End Electronics (GCFE). To cover a large dynamic range of 5 x 10^ 
in each GCFE with commercially available 12-bit ADCs for digitization, the low and high 
energy photodiodes each have their own independent signal chains, the low-energy and the 
high-energy, and each chain operates with two track and hold gains (low and high). This 
arrangement results in four overlapping energy ranges from 2 MeV to 70 GeV overall, as 
shown in Table [51 Range overlap allows cross-calibration of the electronics. Table [5] also 
shows the approximate factors used to convert ADC readout to energy units (MeV). 

A detailed description of the CAL is found elsewhere [8j. 



Name 


Energy 


Gain 


Energy range 


MeV/ADC 


LEX8 


low 


high 


2 MeV to 100 MeV 


0.033 


LEXl 


low 


low 


2 MeV to 1 GeV 


0.30 


HEX8 


high 


high 


30 MeV to 7 GeV 


2.3 


HEXl 


high 


low 


30 MeV to 70 GeV 


20 



Table 5: Four overlapping readout ranges for the CAL and their conversion factors. 

Here we describe the on-orbit measurements of the following: CAL pedestals; crystal 
energy scale, derived from the electronics linearity and crystal light output; crystal light 
asymmetry, which calibrates position measurements along the crystal; and threshold settings. 
These calibrations allow determination of the location and amount of energy deposited in 
each crystal. Processes for estimation of incident photon energy and the resulting overall 
energy resolution are discussed elsewhere [11], [H] . 

5.1 Pedestals 

As discussed in Section 14.11 pedestals are offset voltages for each of the four CAL energy 
ranges that set the "zero point" for the energy scale. We measure pedestals on orbit from 
the periodic triggers issued at 2 Hz by the LAT trigger system during all nominal science 
data acquisitions. Chance coincidence energy deposits result in a small tail to the pedestal 
distribution, but this is suppressed both by the pedestal distribution fitting techniques used 
and by comparison of the various energy ranges. 
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Figure 9: On-orbit pedestal widths (for all channels) for: a) LEX8, b) HEX8, c) LEXl and 
d) HEXl ranges. 



Figure E] shows that typical pedestal widths (RMS) are 0.2 MeV for the LEXl and LEX8 
ranges and 7-10 MeV for the HEXl and HEX8 ranges. The pedestal values are regularly 
monitored on orbit and have been extremely stable since an initial settling period. 

Pedestal measurements made during thermal-vacuum tests in January 2008 indicated a 
channel-dependent linear dependence of pedestal position with temperature, where the drift 
magnitude varies from -3 to +3 ADC units per degree for the LEX8 and HEX8 ranges and 
~10 times smaller for LEXl and HEXl, reflecting their smaller gains. 

When Fermi is in Pointed observing mode, the exposure of the LAT thermal control 
system to the warm Earth changes relative to the Sky Survey observing mode, and the LAT 
detector temperatures change modestly. In the CAL, the temperature changes 1-2 degrees 
within a few hours, leading to a temperature drift of the pedestals of a few ADC units at 
most. Since Log Accept (LAC) thresholds, which determine whether a given crystal readout 
is included in the data stream, are set in absolute ADC unit values (i.e. not relative to the 
pedestal), changes in the pedestal values change the corresponding threshold energies. A 
large enough pedestal change could thus influence data volume. However, only very large 
temperature changes (~ 20 degrees Celsius), which are not anticipated during the mission, 
could lead to significant increases in data volume. 




Figure 10: Characterization of electronics non-linearities. Normalized ADC/DAC versus 
energy for: a) LEXl and b) HEXl ranges. 

5.2 Individual crystal energy scales 

Calibration of the individual crystal energy scales involves the determination of the parame- 
ters of a transfer function between the energy deposited in the crystal and the signal output 
in ADC units. The transfer function consists of the response of the front-end electronics and 
the light emission and collection properties of the scintillation crystals. We use charge injec- 
tion calibrations to describe the electronics and ionization energy depositions from cosmic-ray 
muons or protons and various nuclei to calibrate crystal response. Ionization energy losses 
over a known path length are very predictable and hence make a good calibration tool. 

Charge injection calibrations are used to characterize the non-linear behavior of the 
electronics chain. A pulsed signal of known amplitude (controlled by a charge injection 
calibration DAC) is sent to the preamplifier input of each CAL front-end electronics channel. 
The nonlinearity of this DAC is specified by the manufacturer to be <0.1% and hence 
negligible when compared to the nonlinearity of the front-end electronics. 

For each fixed DAC setting we inject 100 pulses onto each of the electronics channels 
and average the resulting ADC output values. A spline function is used to fit the resulting 
DAC versus ADC curve. The functions (one for each channel), describe both electronics 
gain and non-linearity. Using these functions, signals from each channel can be converted to 
a linear scale. Figure [TO] shows the normalized ADC/DAC versus energy, where deviations 
from one indicate non-linear behavior for measurements made at fixed energy values. The 
largest deviations (~12%) are seen in the LEXl range displayed in Figure [TOk. Figure [TOb 
illustrates the case for the HEXl range with 4% non-linearities. The feature around 2 GeV 
in the HEXl curve corresponds to cross-talk between LEXl, which saturates around 1 GeV, 
and HEXl. For the other ranges (LEX8 and HEX8) deviations are <1%. After applying 
the measured nonlinearity calibration, residual nonlinearity is < 1% of the measured energy, 
resulting in a negligible systematic effect in spectrum determination. 

The crystal response calibration, i.e. the function that relates deposited energy to the 
linearized signal described above, has been performed using different ionizing particles on 
the ground and in orbit. In both cases, ion incident energies are, for the most part, in the 
slow relativistic rise region of the Bethe-Bloch curve, so the predicted energy loss per unit 
path length (^) is only weakly dependent on incident energy. Using simulated incident 
spectra for ground and orbit environments, we determine an expected for each incident 



particle species. We then collect spectra of the these species, correcting each event for path 
length, and compare the peak position in pulse height units to the predicted position in 
energy to yield a calibration. Variations in incident spectrum with orbital position result in 
slight broadening of the energy deposit peaks, but, given integration over multiple orbits, 
the peak most-probable-value is well-determined and usable for calibration. 

On the ground, the low-energy scales were calibrated using sea-level cosmic ray muons 
while high energy ranges were calibrated using muons with the HEXl and HEX8 channels 
set to a special "muon mode" gain setting that increased the gain by a factor of ~10. 

On-orbit, we used a technique we refer to as "proton inter-range calibration". The low- 
energy scales are calibrated using protons. Higher energies are calibrated by using energy 
deposits that meet two criteria: first they must be in the overlap range between LEXl 
and HEX8 and second they must result in a "heavy ion" trigger, the only common trigger 
that produces the required 4 energy range readout rather than the normal single range 
readout. These events are a combination of galactic cosmic ray (GCR) primary carbon 
nuclei, interacting protons and other interacting or ionizing GCRs. From events that meet 
these criteria, we can construct a cross calibration of the low and high energy ranges. In 
both ground and on-orbit cases, the ionization calibration, together with the charge injection 
results yield a usable energy scale that converts ADC units to deposited energy. A week of 
nominal science operations data is sufficient to calibrate the energy scales using relativistic 
protons and heavy ion trigger events in the overlap energy range. 

Protons that are accepted by either the MIP filter (when active) or the diagnostic filter 
(see Table [3]) are required to pass a number of cuts. In order to reject all events that are 
not contained in a single CAL module and to eliminate "corner clipping" events, for which 
path length determination is not sufficiently accurate, we require extrapolated TKR tracks 
to cross the top and the bottom surfaces of a single CAL crystal, and be at least 5 mm away 
from the crystal edges. In addition, we require single TKR track events with > 20 TKR 
hits and a chi-square for the track < 3. To reject low-energy re-entrant albedo protons that 
could broaden and bias the energy distribution, we require the multiple scattering angle, 
calculated by the Kalman filter used for TKR track reconstruction, to be < 0.01. Finally, 
we reject nuclear-interacting protons by selecting events with two or fewer hits in each CAL 
layer and no additional crystals hit in the layer containing the crystal being calibrated. 

In order to determine a calibration peak shape that represents the data well but minimizes 
the number of free parameters, we use a two-step process. First, we produce a spectrum of 
path length corrected signals for all the crystals together. Each peak is fit with a Landau 
distribution convolved with a Gaussian, for which all parameters are left free. From this fit, 
we determine both the Landau and Gaussian widths, which are then fixed. 

In the second step, we fit the Gaussian-convolved Landau function to spectra from each 
crystal separately, allowing peak position and amplitude to vary but using the fixed widths 
determined above. The peak most probable value (MPV) in energy units from the simula- 
tions (10.6 MeV for protons) is divided by the peak MPV in pulse height units determined 
by the fitting process just described to yield the desired calibration quantity for each crystal. 
Figure [11] shows the results, where the most probable value is the peak position. 

Originally, we intended to use "heavy" GCR primary nuclei to calibrate the higher energy 
scales. This would be desirable since 4^ varies as where Z is the atomic number of the 

dx 

incident ion. The major difference between use of GCR heavy nuclei and protons or muons 




Figure 11: Energy deposited in a crystal (pathlength corrected). The position of the proton 
peak is given by the fit. The signal is corrected for electronics non-linearities. 
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Figure 12: Ratio of the energy scales for the on-orbit cahbrations performed in October and 
July 2008: a) low-energy diodes (mean 1.014, sigma 0.009) and b) high-energy diodes (mean 
1.019, Sigma 0.009). 

lies in a phenomenon known as "quenching", in which the crystal light output per unit 
deposited energy (^) is thought to be less for higher Z nuclei than for protons or muons. 
In addition to being Z dependent, this phenomenon also depends on the ion incident energy. 

Since quenching is not well studied for the combination of Z and incident energies relevant 
to on-orbit CAL calibration, we measured the response of the calorimeter CsI(Tl) crystals to 
relativistic nuclei (from carbon to iron) at the GSI facility in 2003 and 2006. The results of 
these studies indicated that for the ion energies examined at GSI, which were rather higher 
than those measured previously for CsI(Tl), ^ was actually higher for measured nuclei with 
6 < Z < 14 than for protons [T^. Due to the lack of a physical model or understanding for 
this "anti-quenching" behavior, we felt that the systematic uncertainties introduced in using 
the heavy ion OCR calibration were unacceptable at this time and have used the proton 
inter-range technique instead. 

In order to study any possible changes in the energy scale calibration between pre and post 
launch measurements and during early on-orbit operations, we calculate, for each crystal, the 
ratio of the energy scale measured after launch to that measured prior to launch. A Gaussian 
fit to this distribution leads to a mean bias, for low-energy diodes, of ~1% and a standard 
deviation which indicates crystal-to-crystal variations of 0.8%. The latter characterizes the 
statistical precision of this calibration procedure. 

For the high-energy diodes, the bias with respect to ground calibration is 5%, which is 
explained by the lack of any high-energy signal to reliably calibrate the ratio between the 
two diodes on the ground, as the muon signal is too small to be visible in the high-energy 
diode with flight gain settings. 

The comparison of on-orbit results from July and October 2008 shows a shift of only 1-2% 
and a spread, crystal-to-crystal, of < 1%. The stability of our calibrations is demonstrated 
in Figure [T2k for the low-energy diode and Figure [T2b for the high-energy diode. 

Despite the fact that they are not suitably well understood for absolute calibrations, we 
do use energy deposits from from GCR heavy nuclei (500 MeV for carbon nuclei and 8 GeV 
for iron) for independent monitoring of the energy scale at high energies. The pathlength- 
corrected spectrum shown in Figure [T3l is obtained by selecting crystal hits in low-multiplicity 
layers, thus rejecting nuclear interactions. Narrow peaks (the carbon peak has a 5% width) 




Figure 13: Energy deposited in all crystals from heavy nuclei collected during 4 days of 
on-orbit operations. Pathlength corrections are applied. 



are easily identified. It is worth noting that the measured energies deposited by cosmic 
rays are consistent with the anti-quenching effect observed in the beam test data acquired 
at GSI. For example, the jp- for carbon is observed to be about 20% higher than expected 
from a Z"^ scaling of the ^ for protons. The count rate in the charge peaks is similar 
to that of the primary galactic cosmic-ray abundance, modified by loss of particles through 
charge- changing interactions above the GAL and by the decreasing efficiency of the on-board 
heavy-ion filter for higher Z nuclei. For example. Figure [T^ shows that, although there is a 
slow systematic shift, the carbon peak position is stable to within 0.5% (for the whole CAL) 
after 2 months of operations. 

5.3 Light asymmetry 

The design of the CAL crystals deliberately "tapers" the light propagation properties of the 
crystal so that we can determine the longitudinal position of an energy deposit by comparison 
of the signal at each crystal end. We refer to the measured quantity as the signal asymmetry, 
defined as the logarithm of the ratio between the signals read out at the two ends of the 
same crystal. 

To calibrate light asymmetry, we select non-interacting heavy nuclei in a similar way to 
protons (see Section [5^ . but with looser TKR track quality. The TKR-determined position 
in a CAL crystal is obtained by extrapolating the TKR track to the center of the crystal 
(half-way through its thickness). Each track enters the crystal through one of the twelve 
evenly spaced bins defined along its length. The distribution of asymmetry signals for each 
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Figure 14: Position of carbon peak for 2 months of on-orbit data. 

bin is collected by computing the asymmetry for each event for the appropriate bin. The 
first and the last bins, near the crystal ends, are not used because of known non-uniformities 
in light collection near the diodes. The average asymmetry is calculated for each of the ten 
central bins along each crystal and is fit with a spline function for purposes of interpolation. 

We use the signals from the two ends of each crystal to obtain the weighted centroid 
of the energy depositions along the crystal. The position resolution along the crystal is 
defined as the RMS of the difference of the position from light asymmetry to that from track 
extrapolation. 

We calibrated the light asymmetry of each crystal on the ground with sea-level cosmic 
muons, and we recalibrated on orbit with GCRs. The calibration constants derived on orbit 
are more precise than those derived on the ground because the GCRs suffer less multiple 
scattering than muons and create larger scintillation signals. The position resolution mea- 
surement for energy depositions from 200 to 900 MeV, for both LEXl, and HEX8 ranges, 
improves from 4 mm to 2 mm and from 13 mm to 9 mm, respectively. Figure [15] shows 
position resolutions with pre- and post-launch constants in the energy range from 200 to 900 
MeV. Results are dominated by the carbon events which peak ~500 MeV. The improvement 
by using flight calibration constants is clearly seen. 

5.4 Trigger thresholds and upper level discriminators 

To reduce the data volume generated by the GAL and the additional dead time that would 
be created by moving large events through the LAT data acquisition system, each GGFE 
has a zero-suppression discriminator with a programmable threshold DAG. This threshold, 
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Figure 15: CAL position resolution from 200 to 900 MeV using ground calibration constants 
a) LEXl and b) HEX8, and using flight calibrations for: c) LEXl and d) HEX8. 
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Figure 16: Signals from one end of a crystal (LEX8 range) and a fit used to determine the 
LAC threshold. 

known as the log-accept (LAC) threshold, is nominally set to 2 MeV, which is approximately 
10 times higher than the average electronic noise. The zero- suppression for an entire crystal 
is performed on the logical OR of the LAC discriminator states at the two ends; thus data 
from both ends of a crystal are included in the CAL data stream if the LAC discriminator 
on either end fires. 

Extensive ground testing with the LAT charge injection system and sea-level cosmic ray 
muons established the functionality, linearity, and energy scale for each LAC discriminator. 
We launched with LAC settings derived from the ground calibrations, but the LAC settings 
are temperature dependent, primarily because the pedestal values are temperature depen- 
dent. Thus we revised the LAC calibration constants on orbit once the LAT had achieved 
its stable operating temperature. For each discriminator, we characterize the relationship 
between DAC setting and LAC value (in MeV) with a linear model that is derived from 
calibration data acquired with the LAC set at two values near the nominal setting. We 
calibrate one end of all crystals at a time using a set of four configurations such that the 
LAC threshold at the crystal end not being calibrated is set to its maximum possible value, 
preventing it from initiating the readout, and the LAC threshold at the end being calibrated 
is set to the test value. This process give a LAC measurement for each channel with a 
statistical precision of 5% (~ 0.1 MeV). Figure [16] shows an example of the signals measured 
at one end of a crystal (LEX8 range) and the fit that determines the LAC threshold. The 
values can be easily converted into energy by using results from Table O We are monitoring 
four GCFEs with out-of-family electronic noise (out of 3072 in the CAL). In February 2009, 
we inhibited one of the channels from participating in zero- suppression decisions because its 



noise level reached 1.5 MeV. This has no impact on the scientific performance of the CAL. 

We monitor the stability of each LAC threshold in all nominal science operations data 
acquisitions. To measure the threshold value, we select events for which we are certain which 
discriminator qualified the crystal for inclusion in the data stream. Because the LAC value 
at each crystal end is within 10% of the setting at the opposite end, we achieve this certainty 
by selecting events where the signal from the two ends differs by more than 10%. 

The CAL provides two fast signals that participate in the formation of LAT trigger, the 
low energy CAL_LO and high energy CAL_HI triggers. The CAL_LO and CAL_HI trigger 
requests are formed as the logical OR of the outputs of the programmable fast low-energy 
and fast high-energy trigger discriminators, respectively FLE and FHE, at each end of a 
crystal. The nominal values for the FLE and FHE thresholds are 100 MeV and 1000 MeV 
per crystal, respectively, as measured at the center of the crystal by each GCFE. 

Extensive testing on the ground with the LAT charge injection system clearly demon- 
strated the functionality and linearity of each FLE and FHE discriminator; however it gave 
only an approximate absolute calibration (i.e. in MeV deposited) of the threshold DACs. 

The FLE and FHE thresholds are calibrated on orbit using background events recorded 
in dedicated-mode with additional information provided by the tower electronics module. 
However, this additional trigger diagnostic information only provides the logical OR of the 
combination of all 12 FLE or all 12 FHE discriminators for each CAL layer-end; thus it does 
not clearly identify which crystal end produced the trigger signal. Distinct procedures for 
FLE and FHE thresholds are necessary to resolve this 12-fold ambiguity. To calibrate the 
FLE discriminators near the nominal 100 MeV setting, we separately enable the trigger for 
each of two groups of six crystals in a layer (the six odd-numbered and the six even-numbered 
crystals) and require that five of the six enabled crystals have a signal below 50 MeV. Most 
showers share their energy between adjacent crystals in a layer, so the separation into two 
sets of six non-adjacent crystals readily resolves that ambiguity. The 50 MeV energy cut 
ensures that only one of the six is the source of the trigger signal. To calibrate the FHE 
discriminators, we allow only the CAL_LO signal to initiate a LAT trigger, and we enable 
the FHE discriminators in two groups of six crystals per layer while we read the state of the 
FHE trigger diagnostic information. 

The efficiency of a discriminator as a function of signal (in MeV) is determined by calcu- 
lating the ratio of the spectrum from events for which the discriminator fires to the spectrum 
for all events. As shown in Figure [T71 the value of the threshold is obtained by fitting a step 
function to this ratio. Having established that the FLE and FHE discriminators are linear 
in threshold DAC setting, we calibrate each discriminator at two settings near nominal value 
and fit the measurements with a linear model. We calibrate FLE at 100 MeV and 150 MeV, 
and we calibrate FHE at 1000 MeV and 1500 MeV. The statistical error in determining the 
threshold values is < 1% for FLE and < 2% for FHE thresholds. Figure [T7I illustrates how 
these efficiencies are obtained for the FLE and FHE thresholds for one GCFE. These data 
are from dedicated calibration runs taken during early operations (July 2008), when thresh- 
olds were set using calibration coefficients derived from ground tests. As it happens, the 
ground calibration gave threshold values for FLE and FHE somewhat higher than intended, 
viz. ~140 MeV and ~1200 MeV, respectively. The ground FLE and FHE calibration relied 
on the charge-injection system, which gives pulse shapes that differ from those produced 
by CsI(Tl) scintillation signals. Since the trigger signal is fast (~ 250 ns) and the energy 
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Figure 17: Efficiency versus energy for the FLE and FHE thresholds. 



measurement depends on the slow shaper (~ 4000 ns), the difference in the shape of the 
pulses is important and creates the 40% and 20% bias we observed. We then adjusted the 
FLE and FHE settings using the calibration constants derived from the on-orbit calibration. 
We continue to monitor the FLE and FHE threshold values with data from the nominal 
science acquisitions. 

The programmable Upper Level Discriminator (ULD) in each GCFE is responsible for 
switching between GAL energy ranges. To select the best range for digitization, three ULDs 
in each GCFE compare the outputs of three ranges (LEX8, LEXl and HEX8) with cor- 
responding ULD threshold (one setting per GCFE). The output of these discriminators is 
analyzed by the range selection logic which selects the best range as the highest range with- 
out an ULD signal fS]. All three ULD thresholds of each crystal end are defined by one DAG 
and are set to ~5% below the saturation level of the ADG. We calibrated the ULD threshold 
DAGs on the ground with charge injection and verified those settings and the linear calibra- 
tion model on orbit during nominal science continuous acquisitions. We measured the ULDs 
on orbit by finding, for each energy range, the highest observed signal for each individual 
crystal end. Most on-orbit data-taking configurations read out only one of the four available 
ranges, namely the one providing largest signal below the ADG saturation, but the nominal 
science data explore all four ranges fully. 

We updated the LAG, FLE, FHE, and ULD threshold settings after their initial on-orbit 
calibrations and later adjusted them to accommodate the small settling drift in pedestal 
values (see Section [5.11 for details). Figure [TS] shows the distribution of all thresholds for 
all channels. The LAG threshold data from August (dashed) and November 2008 (solid) 
are used in Figure [181 where the effect of pedestal evolution since launch is seen as a slight 
broadening of the distribution. This effect is negligible for the other thresholds. The slight 
asymmetry in the distribution of ULD values (Figure [TSH) has absolutely no effect on GAL 
performance; it means only that for a small fraction of channels the range switching will 
happen at a slightly lower energy than for the majority of channels. 




Figure 18: On-orbit measurements of threshold values for November 2008: a) LAC, b) FLE, 
c) FHE and d) ULD. The dashed hne in the LAC histogram shows the effects of pedestal 
drifts seen in August 2008, prior to stabilization. The ULD thresholds are expressed in 
non-pedestal-subtracted ADC units. This makes it easier to judge how thresholds are set. 
The saturation limit corresponds to 4095 ACD units. 



6 TKR calibrations 



The TKR is used to convert the photon to an e"'"/e~ pair and to determine the incoming 
photon direction. It is also the main contributor to the LAT trigger. It consists of sixteen 
modules each composed of a stack of 19 trays. A tray is a stiff, lightweight carbon-composite 
panel with silicon-strip detectors (SSDs) mounted on both sides with strips oriented along 
the same direction. All but the three bottommost trays in each TKR module, contain an 
array of tungsten foils, which matches the active area of each SSD. These foils act as photon 
converters. Depending on their location within the tower (front or back), foils are 3% and 
18% of a radiation length. Each tray is rotated 90° with respect to the one above and the 
one below. Therefore, two consecutive trays are needed to provide orthogonal measurements 
of the x,y coordinates. Details of the TKR design are described elsewhere 

Each side (top or bottom) of the tray consists of 1536 silicon strips read out by twenty four 
64-channel amplifier-discriminator ASICs, GLAST Tracker Front-end Electronics (GTFE), 
which are controlled by two digital readout-controller ASICs, GLAST Tracker Readout Con- 
troller (GTRC). Each channel in the GTFE has a preamplifier, shaping amplifier, and dis- 
criminator similar, although not identical, to the prototype circuits described in The 
amplified detector signals are discriminated by a single threshold per GTFE chip; no other 
measurement of the signal size is made within the GTFE. The TKR electronics is discussed 
in detail elsewhere [7j. 

The trigger information is formed within each GTFE chip from a logical OR of the 64 
channels. Any latched, noisy or inoperable channel can be masked. The OR signal is passed 
to the left or to the right, depending on how the chip is configured, and combined with the 
OR of the neighbor. This procedure continues down the line, until the GTRC receives a 
logical OR of all non-masked channels it controls. This "layer-OR" initiates a one-shot pulse 
of adjustable length in the GTRC, which is sent as a fast trigger signal input for the trigger 
decision. In addition, a counter in the GTRC measures the length of the layer-OR signal, i.e. 
the time-over-threshold (ToT), and buffers the result for inclusion in the event data stream. 
Upon receipt of a signal that acknowledges the trigger decision, each GTFE chip latches the 
status of all 64 channels into one of the four internal event buffers. Another 64-bit mask, 
which is separate from the trigger mask mentioned above, can be used to mask any subset 
of channels from contributing data, as may be necessary in case of noisy channels. 

TKR calibrations include the determination of the noisy channels that form the trigger 
and data masks, of the trigger threshold settings, and ToT calibrations. 

6.1 Noisy channels 

Since noisy channels can increase the false trigger rate, and can affect instrumental dead 
time, ToT measurements and data volume, they are disabled in trigger and/or data masks. 

The trigger mask determines the active channels that can participate in the formation 
of the layer-OR trigger signals. The noise occupancy for each strip is measured using not 
only periodic triggers but all available events, thus increasing statistics by almost 2 orders of 
magnitude. The computation of noise occupancy for non-periodic triggered events excludes 
consecutive (2 or more) layers with at least one hit in each. The occupancy measured by 
this method is consistent with that obtained from periodic triggers. 



Noisy channels produce off-timing trigger signals resulting in a dead time of ~ 1 //s for 
the layer involved. Noisy channels also lead to incorrect ToT measurements if the noise hit 
occurs at the tail of the main pulse. To minimize this effect, we mask any channel with 
occupancy greater than 0.79^. Furthermore, if the occupancy of layer-OR is greater than 
8%, we mask the highest-occupancy channels until we reduce it to this fraction. 

The data mask determines the active channels whose data can be transmitted to the 
ground. Since the offline track reconstruction software is tolerant of high-occupancy channels, 
data mask is driven by the constraint on the data rate given by noisy channels. We mask 
any channel with occupancy greater than 50% since it does not carry any useful information. 
We limit the data size due to noisy channels to be less than 10% of TKR total data size, by 
requiring the average strip occupancy to be less than 5 x 10~^, which corresponds to 44 strip 
hits per event. The typical strip occupancy of ~ 2 — 3 x 10^® is dominated by accidental hits 
due to off-timing cosmic-ray tracks. The strip occupancy due to electronics noise is 10^^ or 
less. We mask highest-occupancy channels until the TKR average occupancy is reduced to 
< 5 X 10"^ 

The number of masked channels for trigger and data purposes was 203 before launch and 
changed to 206 in July 2008, to 220 in August 2008, to 284 in Octobter 2008 and finally to 
316 in January 2009. These additional 113 channels are distributed across seven SSDs, while 
60 of these channels are concentrated in a region of one SSD. The total number of disabled 
channels corresponds to only 0.04% of the total number of TKR channels. 

6.2 Trigger and data latching thresholds 

In order to minimize the noise occupancy while maximizing the hit efficiency, the nominal 
threshold level is set to 1.4 fC (~0.28 MIP). The threshold DAC value for each GTFE is 
calibrated using charge injection. The charge injection DAC is set to the value corresponding 
to 1.4 fC and the threshold DAC is scanned (see Sectior i^74]) . The best threshold for each 
channel is determined by a fit to the occupancy versus threshold using the error function 
(integral of a Gaussian). The average threshold for each GTFE is obtained by calculating 
the mean value of the threshold DAC values after removing all dead and masked channels, 
and 5% of the channels with the largest and the smallest values. 

The calibration of the threshold DAC performed after launch yields identical DAC values 
to those before launch for 86% of GTFEs. Only 0.04% of the total number of GFTEs (6 out 
of 14 000) exhibits a difference of more than one DAC value. This was already known from 
pre-launch measurements and it corresponded to additional noise in the system. 

The GTFE data is latched following a trigger acknowledge signal, which is 0.8 /is later 
than the typical peak of the TKR pulse shape. This value was determined using an external 
trigger during pre-launch testfl Due to the delayed data latch timing, the effective threshold 
for the data is different from the trigger threshold. Once the threshold DAC value is deter- 
mined for all GTFEs, we measure the effective thresholds at the time of the data capture 
by scanning charge injection calibration DAC values. The best threshold for each channel 

^0.7% corresponds to the value in which the loss in efficiency due to dead time is comparable to that fi'om 
masking noisy channels. 

^Although this implies that the TKR delay should be negative, the smallest allowed value is zero. This 
has no effect on the TKR performance. 



is determined by a fit to the occupancy versus threshold data using the error function. The 
ratio of data latching threshold measurements before and after launch yields a 2% shift in 
RMS and 0.5% shift in the mean value, which implies no significant changes from values 
measured prior to launch. Therefore, trigger thresholds were assumed not to have changed 
after launch and were not recalibrated. 

Figure [12] shows the trigger thresholds obtained prior to launch and the data latching 
thresholds measured on-orbit. Because of the delay in latching the data, the threshold for the 
data capture is slightly higher and has a broader peak than that of the trigger threshold. The 
RMS dispersion is ~5% and ~12% for the trigger and data latching thresholds, respectively. 
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Figure 19: Effective threshold for data capture (solid) and trigger threshold (dashed). 



6.3 Time-over-threshold conversion parameters 

To determine the conversion parameters from ToT (ns) to charge deposit (fC), we measure 
for each channel, the ToT value for several settings of the charge injection calibration DAC. 
Since by definition the ToT values cannot be negative, the ToT measurements near threshold 
are biased toward positive values and result in slightly biased conversion parameters. Fig- 
ure EO shows the amplitude of charge injected versus ToT values, where the values for charge 
injection have not been corrected by the MIP scale calibration described in Section 16. 4[ In 
the fit shown in Figure [20l the ToT is described as a second order polynomial of injected 
charge, whose offset corresponds to the threshold value. To avoid biases near the threshold, 
the fit uses a fixed value for the intercept that corresponds to the calibrated threshold. The 



statistical error in the fit of order 8% and is estimated by comparing two measurements of 
the same curve done prior to launch. This calibration is important since the ToT gain may 
vary by as much as a factor of three within a GTFE. We have not performed this calibra- 
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Figure 20: Amplitude of charge injected versus ToT values. The curve corresponds to a 
second order polynomial fit to the data. 

tion after launch, since other TKR results (see Section 16.20 indicated little change in the 
response of TKR pulses and this calibration requires more than 10 hours of data-taking in 
dedicated-mode. We expect to recalibrate these values annually. 

6.4 MIP scale calibration 

The absolute calibration of the charge injection DAG was performed prior to launch using 
the charge deposited by surface cosmic rays. A correction factor was defined by the ratio 
of MIP peaks between data and Monte Carlo simulations. The distribution of this ratio 
shown in Figure [2T] exhibits an RMS dispersion of ~9%. The MIP deposited in each channel 
was calibrated using these correction factors. On-orbit, we select single-track events, close 
to normal incidence (cos^ < —0.85) and with CAL energy consistent with that of a MIP. 
We require the RMS values of TKR hit positions with respect to the reconstructed track 
position to be consistent with the resolution needed to reject low-energy tracks. To avoid 
confusion with charge sharing between adjacent strips, we only consider layers with single 
hit strips. The charge deposited by a particle is corrected by taking into account its path 
length in the silicon. The measured ToT of the hit associated with the track is converted 
to charge using the conversion parameters. Following that, we fit the data for each GTFE 
with a Landau distribution convolved with a Gaussian. 




Scale correction factor 

Figure 21: Measured MIP peak divided by simulated MIP peak obtained prior to launch. 
This charge scale correction factor is obtained for all GTFEs. 



More than 10 million MIP tracks are required to accumulate sufficient entries for all GT- 
FEs, which takes 5 days during nominal science operations. Figure [22] shows the MIP charge 
deposit distribution for all channels before (dotted histogram) and after (solid histogram) 
each GTFE is calibrated. The dispersion correction factor of 9% is included in the calibrated 
results. 

7 Determination of SAA polygon 

The orbit of Fermi intersects the Earth's inner radiation belt in a region which is known as 
the South Atlantic Anomaly (SAA). This region features geomagnetically trapped protons 
with energies up to hundreds of MeV and electrons with energies up to tens of MeV. The 
flux of protons and electrons in the LAT energy range reach levels which are several orders 
of magnitudes above those of primary cosmic rays. This extreme particle flux imposes 
constraints on LAT operations. The TKR electronics saturate due to the increase in the 
charge deposited per live time, leading to large dead time fractions, thereby hampering 
scientific observations. The continuous infiux of particles generates high currents in the 
ACD photomultiplier tubes (PMT), thus exceeding safe operating limits, which leads to 
slow deterioration. Therefore, during SAA passages, triggering, recording and transmission 
of science data are stopped and the bias voltages of the PMTs are lowered from 900V to 
~400V. Only LAT housekeeping data are recorded and transmitted to the ground. 

The position along the orbit defined by the GPS receiver aboard the Fermi spacecraft 
determines the transition between nominal science operations and the SAA transit mode. 
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Figure 22: MIP signal (pathlength corrected) before (dotted) and after (solid) calibrations. 



The latitude and longitude of the Fermi position are compared to the bounds of a polygon 
defined by 12 latitude-longitude vertices stored in the spacecraft memory. As the spacecraft 
position crosses this polygonal boundary it triggers the SAA transit mode. To avoid multiple 
entries and exits during a single orbit, a convex polygon is used to define the SAA region. 

We chose a conservative definition for this initial SAA boundary, with the expectation 
that we would update the boundary based on particle rate measurements made with the LAT 
once it was on orbit. The first version of the polygon, or SAA boundary, was determined be- 
fore launch based on models of the Earth radiation belts and data from other spacecraft. The 
inner radiation belt was modeled using trapped radiation models: AP-8 [18] and PSB97 [19] 
in conjunction with the current version of the International Geomagnetic Reference Field 
(IGRF-10) [20]. The 12-edge SAA boundary polygon was calculated from these models 
based on the contour in latitude and longitude, where the E > 20 MeV trapped proton flux 
reached 1 cm~^s~^. For regions where two models predictions disagreed, we chose the larger 
flux. The smallest convex polygon circumscribing this contour was selected and padded by 
a 4° margin. Figure [23k shows the trapped proton flux profiles above 20 MeV predicted by 
trapped radiation models versus geographic latitude and longitude. Figure. [23b shows top 
scintillator count rates in the altitude range between 532 km and 575 km, reported by the 
PAMELA experiment |21], which is similar to that of the Fermi orbit. The red polygon 
shows the locations of the SAA boundary edges determined before launch and used during 
the initial phase of the Fermi LAT mission. This definition of the SAA boundary resulted 
in a loss of observation time of about 17%. 

After launch, diagnostic data of the LAT were used to refine the size of the polygon. 
Even though science triggers are disabled during SAA passages, fast trigger signals remain 




Figure 23: Trapped radiation models versus geographic latitude and longitude: a) trapped 
proton flux above 20 MeV predicted by the AP-8 [TH| (black contour lines) and PSB97 IL9\ 
(color contours), b) top scintillator count rates reported by the PAMELA experiment [2T]. In 
both of the plots, the red polygon shows the locations of the SAA boundary edges determined 
before launch and used during the initial phase of the LAT operations. 



operational. Special TKR and ACD counters can sample the rate of fast trigger signals to 
determine position-dependent rates of the LAT along the orbit. Figure [21] shows the rates 
recorded in the TKR counters versus spacecraft position. A rate increase is visible at the 
edges of the SAA before the TKR electronics saturates and suppress fast trigger signals, 
thus bringing the count rates to zero. 

To define the SAA boundary using these data one has to account for particle rates and 
rate variations from primary and secondary cosmic rays, where both depend on the local 
geomagnetic cutoff rigidity at the Fermi location. Therefore, we used data from the region 
outside the pre-launch SAA boundary to determine these rates and set an upper limit on 
the expected number of cosmic-ray counts per second as a function of the local geomagnetic 
rigidity cutoff. The optimized SAA boundary polygon is calculated including the points 
that exceeded this limit. A padding of 1° is applied to account for the limited resolution and 
sensitivity of the measurement. 

Figure [21] shows the average rate of TKR counters obtained during 26.6 days of LAT 
nominal science operations versus geographic latitude and longitude. Superimposed are the 
pre-launch SAA boundary (red) used during the initial phase of the mission, and the refined 
polygon (yellow) uploaded to spacecraft memory. The updated polygon reduced the loss in 
observation time to approximately 13% of the total on-orbit time. This polygon has been 
the default for the LAT operations since July 28, 2008. A cross-check during nominal science 
operations is performed with the ACD trigger signal counters. These are more sensitive to 
the low-energy component (E ^ 10-60 MeV) of the trapped particle flux than those from 
the TKR. There is no significant increase in the rate of ACD fast trigger signals as Fermi 
approaches the SAA boundary, thus validating the optimized polygon. 

Since the SAA moves at a rate of a few tenths of a degree per year and its size and 
particle fluxes vary with the solar cycle, we expect annual updates to the SAA boundary. 
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Figure 24: Average rate of TKR counters obtained during 26.6 days of LAT nominal science 
operations versus geographic latitude and longitude. Superimposed are the pre-launch SAA 
boundary (red) used during the initial phase of the mission, and the updated SAA boundary 
(yellow) derived from measurements of the TKR counter data (see text) 

8 Live time 

The live time is accumulated taking into account the variety of dead time effects. Science 
data taking is disabled during SAA passages (see Section [7j). Instrumental dead time, during 
event latching and readout, is about 8% on average outside the SAA, although this fraction 
depends on the trigger rate (primarily background). The dependence on the geomagnetic 
latitude of Fermi is secondary. Other losses are caused by failures in transmission and ground 
processing and dedicated-mode calibrations. 

Accurate accounting for live time is essential for obtaining calibrated fluxes and spectra 
for astrophysical sources of gamma rays. The live time relates the effective collecting area 
of the LAT to the overall exposure. Owing to the very large field of view of the LAT 
(>2 sr) and the relatively slow scanning rate (~4 deg min~^) the accumulated live time 
typically is needed only coarsely (~30 s intervals) for accurate exposures to be calculated. 
For very bright transient sources, finer accounting for live time is used, owing to the high 
and variable rates of triggers. For example, a bright GRB in the field of view of the LAT, 
such as GRB08019C [22], can induce a dead time fraction of about 16% during the impulsive 
phase of the burst. 

The dead time for event latching and readout is tracked on-board the LAT every 50 ns 
using the 20 MHz LAT system clock. The stability of the 20 MHz clock is closely tracked 
with a 1 pulse-per-second (PPS) signal from the GPS system of the spacecraft. 

During nominal science operations, the instrumental dead time is dominated by the fixed 
time to read out the event timing and trigger information, which imposes a minimum dead 
time of 26.5 //s per event. Every time the LAT triggers, further data taking is disabled until 
the data from the event are read out. After the end of the time coincidence window the 



latency of the trigger system is 100 ns before the new time coincidence window is available, 
even if the previous event was not read out. Single front-end electronic channels can be dead 
for several microseconds while the pulse is above threshold. In special cases, such as when 
the periodic trigger is enabled (2 Hz), the entire CAL is read out with no zero suppression 
and in this case, the dead time can be as large as 620 /is. During nominal science operations 
one of the trigger combinations is configured to read the CAL with all four energy ranges 
and zero suppression for which the dead time is about 65 fis. 

Losses of data in transmission can effectively remove events from the data stream both 
directly and indirectly. The indirect effects relate to how the events are assembled in the 
telemetry stream. In each packet, the event times are encoded as times relative to the 
start of the packet. In addition, the 20 MHz counters mentioned above roll over every 1.6 
s and the times of the roll overs are also encoded in the packets. Loss of certain parts of 
a packet can cause times and live time accumulations to be lost for up to ~200 events. 
At the 450 Hz nominal event rate in telemetry, these losses can be as large as 0.4 s. The 
mission specification for data loss is less than 2%, and in practice loss due to transmission 
errors has been much less than this, in part because on-board the LAT data are retained 
for approximately 24 hours and retransmissions can be requested. At the ground station, 
Fermi telemetry is buffered for 1 week and losses in ground transmission generally can be 
recovered. The overall average loss has been much less than 1% to date. 

Losses in ground processing of event data are also rare and are significantly less than 1%. 

9 Overall timing accuracy 

Recording accurate arrival times of LAT photons is essential for studies of gamma-ray bursts 
and pulsar timing. Absolute timing tests were performed during pre- and post-launch activ- 
ities. A discussion of LAT pulsar timing can be found elsewhere [23] . 

During pre-launch tests we recorded cosmic rays to measure the time difference between 
two GPS systems. As shown in Figure [25] a pair of scintillator tiles provided a reference 
for the LAT timestamps. The coincidence signal from these tiles triggered a VME-based 
GPS time system previously used by the ground gamma-ray telescope CELESTE. Its abso- 
lute time accuracy was previously demonstrated by measuring the Crab optical pulsar [21]. 
Reconstructed muon tracks traversing the LAT detector were extrapolated to their impact 
point on the laboratory floor and their timestamps were measured with respect to the GPS 
of the Fermi satellite. If a muon passed through the pair of scintillators placed next to Fermi, 
a GPS timestamp from a standalone VME data acquisition system was also recorded. Fig- 
ure [26] shows that the LAT timestamps agreed with the reference GPS to within 0.3 /is. 

GPS receivers use the arrival times of reference signals from other GPS satellites to 
calculate their time and position and transmit that information to processors on Fermi. 
This is accompanied by an electronic "Pulse Per Second" (PPS) at the moment of validity 
of the timestamp word. The processors, using precision oscillators on-board the spacecraft, 
maintain the PPS accuracy in the case of occasional short losses of GPS signal reception. 

Bright gamma-ray pulsars were used to verify that the integer seconds of absolute time 
from the GPS receiver conform to Coordinated Universal Time (UTC), since an integer offset 
in the Fermi clocks would make a large shift in observed gamma-ray phase, different for each 
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Figure 25: Diagram of the muon scintillator telescope placed next to the Fermi satellite 
during pre-launch tests. 



pulsar. The rotational phase of the gamma-rays peaks of the Vela pulsar relative to the radio 
peak agree with that measured by previous experiments [25j. The first gamma-ray peak of 
the Crab pulsar is aligned with the radio reference, since Crab radio and gamma-ray beams 
appear to come from the same part of the neutron star magnetosphere [23] . 

The fractional part of LAT event timestamps is obtained from the counts of a 20 MHz 
oscillator recorded by scalers latched at the reception of a GPS PPS signal and at the 
reception of an event trigger. The behavior of the oscillator was extensively characterized 
during the ground tests, and its frequency is recalibrated each second using the PPS-latched 
scaler values. On-orbit telemetry monitoring shows that the internal spacecraft timing signals 
behave as before launch, from which we conclude that LAT timestamps are still well within 
1 fis of the GPS times used by the spacecraft. GPS times are maintained within 20 ns (1 
Sigma) of UTC [2S]. 

At present the best observational validation of the on-orbit clock performance comes 
from the pulsars PSR J0030+0451 and PSR J1939+2134. The peak width of < 120 fis 
reported for PSR J0030-I-0451 in [27] demonstrates the stability of the LAT event times over 
six months of data-taking, but not their absolute accuracy. However, the 1.56 ms pulsar 
PSR J1939-I-2134 appears to have a gamma peak aligned with the radio peak to better than 
1/20 of a rotation of the neutron star, that is, 80 ^s [28j. This translates to an absolute 
time accuracy if one assumes that both the gamma and radio emission come from the same 
region in the pulsar's magnetosphere. 
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Figure 26: Results from pre-launch tests: a) differences in the times recorded witli tlie two 
acquisition systems, versus elapsed time; b) histogram of time differences indicating mean 
and RMS values around 0.3 fis. 



10 Internal and spacecraft boresight alignments 

The accuracy of the directions of reconstructed events dependes on the knowledge of the 
exact position of each element of the TKR, i.e individual SSD or planes. Misalignments of 
any element can broaden the instrument response function, thus lowering the LAT sensitivity. 
We perform the following alignment procedures: 

1. intra-tower alignment to determine the position and orientation of each element of a 
single tower with respect to an ideal coordinate system; 

2. inter-tower alignment to determine the position and orientation of each tower with 
respect to an ideal coordinate system; 

3. spacecraft alignment to determine the rotation of the LAT with respect to the Fermi 
on-board guidance, navigation and control system. 

10.1 Intra-tower alignment 

As described in Section [6], each TKR tower consists of 36 silicon planes each instrumented 
with sixteen silicon microstrip detectors. All strips in a plane are oriented along the same 
direction. The intra-tower procedure aligns the planes horizontally (along the measured 
coordinate) and vertically (perpendicular to the silicon plane), and determines all rotations 
of the planes within the TKR tower. 

Due to the procedures for construction and assembly of towers, we expect only minor 
deviations from the original positions. This simplifies the formalism used for the alignment, 
and we assume first-order approximations to deviations. 

The trajectory of a particle traversing the TKR is characterized by a reconstructed track, 
which consists of a list of associated silicon strip hits joined by a straight line. As the particle 



Figure 27: The diagram shows five sihcon layers (xO to x4) and the ideal x2 layer, crossed by 
two tracks. Measured hits are denoted as filled diamonds, extrapolated hits (assuming the 
ideal position of the x2 plane) as open diamonds. The residual is composed of two contri- 
butions: Ax due to horizontal misalignment, and A2; ■ tan(0) due to vertical misalignment. 

encounters material in its path its direction is affected by multiple scattering. The effect 
is more pronounced in the presence of dense materials such as the tungsten converter foils. 
Thus, a fit to an otherwise straight track results in deviations of the real hit positions from an 
ideal straight track referred to as residuals. The diagram in Figure [271 shows two tracks in the 
xz plane and illustrates how residuals relate to measured and ideal positions. A distribution 
of residuals should be centered at 0, and any deviation from zero is an indication that the 
element studied is not at its assumed position. Figure [28] illustrates how straight tracks 
are used to determine the horizontal and vertical misalignments. It shows the hit residuals 
versus tan(^) (i.e. inverse of the slope) of the track, for one silicon plane (arbitrarily chosen). 
A straight line fit through these points produces an offset and a slope, which correspond 
to deviations of the real position from the assumed one. Results from the fit determine 
horizontal (Ax) and vertical ( Az) deviations, corresponding to shifts along and perpendicular 
to the strip coordinates, respectively. 

Intra- and inter-tower TKR alignments are iterative processes that use reconstructed 
events recorded during nominal science operations with no requirement on any on-board fil- 
ter. As a result, each TKR element is subject to a wide range of event types, thus minimizing 
the possibility of selection bias. 

Here we briefiy describe a method common to both intra- and inter-tower alignments. 
The event reconstruction associates hits to tracks and classifies them according to their 
track length and straightness. The internal TKR alignment procedures benefits from track 
reconstruction by using only the first (best) track in the event, i.e. all other tracks are 
discarded. The list of the hits is then input to the alignment algorithm. 

The alignment algorithms rely on two important pieces: the raw position information, 
i.e. the plane and strip containing a hit, and the smoot/ied position, which results from track 
finding and optimization through the Kalman filter algorithms [I]. Since the raw position 
cannot be used to determine the location of the hit along the strip, the algorithm uses the 




Figure 28: Residual versus tan(^) of the track, for a silicon plane arbitrarily chosen. A fit 
with a straight line provides an offset and a slope, which correspond to deviations of the real 
position to the assumed one. 

corresponding smoothed position, instead. Finally the algorithm fits the hit positions with a 
straight line, computes the of the fit, and evaluates the residuals, and the angle 6 of the 
track with respect to the tower axis. 

Correlations similar to those depicted in Figure [28] are used to determine all alignment 
parameters listed in Table El Note that the angle 6 is always measured in the coordinate of 
the silicon plane (e.g. 6x for x-planes). The sign of the correlation for rotz is different for 
the X and the y-planes. There are six parameters to determine, but only four correlations. 
Each correlation gives two parameters: the horizontal position, either Ax (for x planes) or 
Ay (for y planes), and Az or any of the three rotations. For each TKR element the iterative 



Correlation of track parameters 


Alignment Parameters 


residual versus tan(6') 


Ah + Az 


residual versus x ■ tan(^) 


Ah + rotx 


residual versus —y ■ tan(^) 


Ah + Yoty 


residual versus =F position along the strip 


Ah + Totz 



Table 6: Correlations of track parameters, and resulting alignment parameters. Ah is the 
measured horizontal coordinate of a plane, i.e. Ax for x-planes and Ay for y-planes. See 
text for more details. 

algorithm determines five parameters that are measured relative to other elements whose 
positions are also unknown. The knowledge of the precise positions and rotations increases 
with successive iterations, and each parameter converges to its correct value. The decision 



for convergence comes from comparing the largest deviation of any positional (rotational) 
parameter between two successive iterations and by requiring the maximum deviations to be 
0.01 fim (1 /irad). Convergence typically occurs after 100 iterations. Currently, the smallest 
elements considered in intra-tower alignment are single silicon planes (top and bottom of 
trays as defined in Section [6]). In the future, we will extend the method to align SSD's. 

Table [7] shows the results from the on-orbit intra-tower alignment of 576 TKR planes. 
Offsets and slopes (as in Figure 1251) are entered in a histogram and the standard deviation 
(a) of each distribution is shown in Table [71 Since planes are aligned with respect to an 
ideal frame, by construction, the sum of all shifts (for all parameters) average to the mean 
value of 0. Results clearly demonstrate the quality of the assembly of individual towers since 
standard deviations are within ±61 and ± 220 /irad. 



Parameter 


Standard deviation (a) 


Ax 


± 43 fim 


Ay 


± 59 fiia 


Az 


± 61 fim 


rota; 


± 220 //rad 


roty 


± 220 /irad 


rot^ 


± 210 /irad 



Table 7: On-orbit intra-tower alignment constants, for 576 planes (averaged for all towers). 



10.2 Inter-tower alignment 

The inter-tower alignment procedure aligns towers spatially and determines their rotations 
with respect to the LAT reference frame, and is performed only after the intra-tower align- 
ment constants have been determined. While the latter uses the residuals of each hit to 
determine the alignment constants, the former relies on track segments of events that cross 
tower boundaries. We determine the orientation of each tower by evaluating the angles be- 
tween pairs of track segments. The rotation in one axis is given by the scalar product of the 
track segments, projected on the plane spanned by the two other axes, and averaged over 
all selected events: 

cosM,) = ^5^^^, (1) 

where t and ? are the track segments in vectorial form and i corresponds to x, y, or z. For 
convenience, Eq.([l]) can be rewritten in terms of track slopes cot{9), yielding 

After fixing rotations, each pair of track segments is projected onto a reference plane 
between both towers, as depicted in Figure [291 Each intersection defines a position on this 
plane. For perfect alignment, these positions are identical. The differences between both 
positions are analyzed as in Section 110.11 i.e. by searching for correlations between the 
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reference plane 

Figure 29: Horizontal and vertical alignment for a pair of towers, the real positions are 
drawn with solid lines, the assumed position of tower B with dashed ones. The black line 
between the towers marks the reference plane. The real particle track is shown with solid 
line, the extrapolated track based on the assumed position dashed. 

position differences and cot(^) of the track. In the example of Figure [29l the residual has a 
Az and a Ax- cot(^) contribution due to vertical and horizontal misalignments. 

Table M summarizes the results from the on-orbit inter-tower alignment of 16 TKR tow- 
ers. Contrary to the intra-tower alignment, only a small number of values (16 towers instead 
of 576 planes) are available to determine the a of the distribution for each of the alignment 
parameters. Because of that. Table [8] has an additional column that shows the largest and 
smallest value for each parameter. Results clearly demonstrate the quality of the assembly of 
towers in the LAT since standard deviations are within ± 120 fim and ± 260 /irad. Alignment 



Parameter 


Values (mix,max) 


Standard deviation [a) 


Ax 


(-250,+190) fim 


± 119 ^m 


Ay 


(-90,+130) fim 


lb 68 /im 


Az 


(-150,+150) fim 


± 87 /im 


rotx 


(-450,+400) ^rad 


± 260 /xrad 


Toty 


(-360,+480) //rad 


± 250 niad 


TOtZ 


(-360,+360) firad 


± 230 niad 



Table 8: On-orbit inter-tower alignment constants for 16 towers. 



procedures were validated using a simulated dataset of cosmic protons, generated with ran- 
domly misaligned geometry for the LAT. Figure [30] shows the results when the alignment 
procedure is applied to data acquired before and after launch. Here the standard deviations 
(cta) are calculated from a histogram of differences between ground and on-orbit measure- 
ments. Figure [5Uh shows the correlation for intra-tower positional alignment constants, while 





ground position (|am) position differences (urn) 

Figure 30: On-orbit data: a) correlation between ground and on-orbit inter-tower positional 
alignment constants, b) the histogram of the differences (Aa), where x, y and z correspond 
to squares, diamonds and triangles, respectively. 

Figure [30b displays the histogram of these residual differences. Table |9] compares unbiased 
datasets obtained on ground (~8 million surface muons) with orbit data (~20 million events, 
mostly cosmic protons) and lists the average differences for both intra- and inter- alignment. 



Parameter 


(Ta (intra-tower) 


(Ta (inter-tower) 


Ax 


2.4 fim 


30 /im 


Ay 


2.5 fim 


32 /im 


Az 


8.8 fim 


17 /im 


rotx 


51 /irad 


53 /xrad 


Toty 


53 yurad 


42 /xrad 


TOtZ 


36 /irad 


24 /xrad 



Table 9: Standard deviations (cta) of the differences between ground and on-orbit alignments. 
Inter-tower values are averaged over all three positions and rotations. 

for all parameters. Results from Figure [30] and Table [9] demonstrate that misalignments due 
to launch or temperature variations on-orbit are small and limited to within ± 35 /xm and 
± 55 /irad. 

Finally, Figure [31] shows all correlations between the rotation angles obtained for each 
tower using on-orbit data acquired by the end of October 2008. Tower numbers are labeled 
from to 15. Squares denote the correlation between rotx (x-axis) and roty (left y-axis label), 
and diamonds between rotx and rotz (right y-axis label). The spread in the correlation 
between Totx and roty translates directly into a pointing inaccuracy on the celestial sphere, 
e.g. a misalignment of 0.3 mrad in rotx or roty for one tower would cause a source position 
error of the same size. The absence of clustering in this distribution indicates that there are 
no obvious systematic effects. 
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Figure 31: Correlation between the different rotation angles of each tower (the number of 
each tower is denoted in the diagram). Squares denote the correlation between rotx (x-axis) 
and roty (left y-axis label), and diamonds between rotx and rot^; (right y-sods label). 



The alignment of the TKR has been checked twice before launch (since October 2006) 
and four times after launch (since June 2008). No significant changes of the positions and 
rotations were observed. We expect updates to be made annually. 



10.3 LAT alignment with respect to the spacecraft 

The LAT reconstructs the direction of each photon in the LAT reference system. To convert 
this to celestial coordinates, we require the celestial orientation of the LAT. The spacecraft 
orientation is provided by the Guidance, Navigation and Control system, or star-tracker. 
This system involves optical telescopes on an optical bench, and a star field pattern recogni- 
tion system, which was calibrated before launch. The LAT system is nominally the same as 
the spacecraft system, but small deviations are expected due to uncertainty in the ground 



alignment process, thermal variations, launch vibrations, relaxation in Og, or the inter-tower 
alignment process. Measuring these deviations, by comparing the LAT view of the sky with 
that seen by the star-tracker is called "boresight alignment" . 

The boresight alignment is determined via likelihood maximization for gamma rays near 
bright, identified celestial point sources of gamma-ray emission. The likelihood analysis is 
"binned" in the sense that the gamma rays are not considered individually but binned into 
maps in celestial coordinates. Separate maps are used for different energy ranges because the 
point-spread function (PSF) depends strongly on energy. The likelihood is defined in terms 
of the Point Spread Function (PSF), the background (sum of celestial sources and residual 
charged particle bacgrounds, estimated independently for each source), and the boresight 
alignment angles. We measure the boresight alignment by maximizing the total binned 
likelihood of the point source model (background is included) with respect to three angles 
characterizing an arbitrary rotation. We chose the convention of performing consecutive 
rotations about the x, y, and z axes. These are defined by the spacecraft: z is along the 
symmetry axis (pointing direction), and y is along the solar panels. Since the angles are 
small, less than a degree, the order of the rotations is irrelevant. Rotations about the x and 
y axes thus correspond to angular deviations of the same size, while the rotation about the 
z axis affects only tracks that traverse the detector at off-axis angles. The procedure for 
calculating the likelihood transforms both the source position (i.e. celestial coordinates) and 
the photon position into spacecraft coordinates, rotates the photon's position, and determines 
the likelihood that the photon is consistent with the source location. Details of the likelihood 
method are discussed elsewhere [2S]. 

We first define a set of reference bright point sources using the likelihood itself to select 
the most significant. We select photons within the 99% PSF containment of the sources. 
Since the angular containment is energy dependent, photons are divided into energy bands, 
using an average PSF for each band to define the expected deviation distribution. The 
PSF is non-Gaussian, and has been found to be well modeled by the following power-law 
expression [TT] . 



where Ei is the energy bin of the photon, 5 is the angular separation between recon- 
structed and true direction, and 7(i?i) and a{Ei) are energy dependent parameters. The 
parameter j{Ei) determines the tails of the distribution and at low energues has a value of 
~2.25, which decreases at high energies, creating longer PSF tails. The parameter a{Ei) has 
a power-law dependence at lower energies and reaches an asymptotic value at high energies 
determined by the silicon detector readout pitch. We compute 6 for all spatial bins and 
weight the likelihood by the number of counts per bin. The parameters of the PSF were 
estimated by extensive Monte Carlo simulations, but have been since refined with on-orbit 
data. We select photons from energy bands above 500 MeV, which contain little or no back- 
ground. We calculate a likelihood for each source by modeling it as point source in a uniform 
background. Treating the remaining background on the angular scale of the PSF as uniform 
is a good approximation, and would be significant only in the tails where the contribution 
to the likelihood is small. 




(3) 



The bright point sources for daily ahgnment check are selected based on the Test Statistic 
(TS), which is twice the difference of the log likelihood of a point source in a uniform 
background and the log likelihood of a purely uniform background for a particular energy 
band. The total TS is then the sum of the TS's for the energy bands above 500 MeV. Since the 
background fraction is already small, the assumption of uniformity is a good approximation 
to the distribution of the background events, and we select all sources with a total TS greater 
than 25. For the example shown in Figure [32], 6 sources match this criterion, and these are 
used as seeds for the alignment procedure. The largest contribution comes from the Vela 
pulsar, which is dominated by photons between 1 GeV and its 10 GeV cutoff. 

The likelihood is maximized with respect to the background for each energy band, since 
the shape of the PSF depends on energy. The total likelihood for a point source is then the 
product of the likelihood in each energy band. 

Figure [32] shows typical projected error ellipses resulting from the likelihood fit. The 
crosses correspond to the ideal location and the dashed and solid contour lines represent 1 
and 2 sigma contours derived from the log likelihood, respectively. A typical day in orbit 
yields a precision of 1.5 arcmin in x and y, and 2.4 arcmin in z, when using 174 photons 
collected from six astronomical sources. We expect, after the analysis of the first year data 
to reach the required value of 4 arcsec for the boresight alignment residuals. 

As shown in Figure [33] we make independent measurements of the boresight alignment 
for each week, to monitor the stability. The figure shows the cumulative mean (dash-dotted 
line) for the rotation angles about the x, y and z axes for a period of 3 months. The 
cumulative mean is just the average of the parameters weighted by the errors. We use 
the same sources and data to optimize the combined likelihood with respect to the PSF 
parameters themselves, and feed this back into the analysis. Weekly measurements indicate 
a stability of 0.3 arcmin, constant with the statistical errors. 

We repeated the likelihood fits for a period of about 5 months. The results for the rotation 
angles are shown in Figure [311 We clearly see the improvement with more statistics. As 
expected, rotations about x and y axis lead to similar results (6 arcsec), while for rotations 
around z we obtain values approaching 8 arcsec. These values should be added in quadrature 
to other contributions to obtain the systematic uncertainties for localizing point sources. 

11 Conclusions 

We have discussed the on-orbit calibrations for the Fermi LAT, which include synchroniza- 
tion of trigger signals, optimization of delays for latching data, determination of detector 
thresholds, gains and responses, evaluation of the perimeter of the South Atlantic Anomaly 
(SAA), measurements of live time and of absolute time and internal and spacecraft boresight 
alignments. The results summarized in Table [TU] were obtained using known astrophysical 
sources, galactic cosmic rays, and charge injection. There were only minor changes to cal- 
ibration constants since launch and these quantities have been stable during the first eight 
months of operations. 

The LAT with almost a million channels is a remarkably stable instrument and is ex- 
pected to operate during the next few years with dedicated calibration runs of about three 
hours every three months. The frequency with which calibration updates are expected to 
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Figure 32: Projection of the likelihood surfaces into each of the planes a) yx, b) zx and c) 
zy for one day of data (August 3, 2008). Crosses correspond to the ideal location (perfect 
alignment), and contours are 68% (dashed) and 95% (sohd) containment. 
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Figure 33: Results of the likelihood fit for weekly measurements over a period of 3 months 
for the rotation angles about the a) x, b) y and c) z axes. The cumulative mean is displayed 
as dash-dotted hue, while the dashed line shows the reference at zero {y axis scales are 
different). The start date for these measurements on Mission week 10 corresponds to August 
10, 2008. 




Figure 34: Data correspond to results of the likelihood fit for measurements over a period 
of 5 months for the rotation angles about the a) x, b) y and c) z axes. 



occur varies from three to twelve months, but these do not necessarily impact the LAT per- 
formance in any significant way. As a consequence, changes to high level datasets due to 
these calibrations will be infrequent. The results reported here have been used to calibrate 
the LAT datasets to be publicly released in August 2009. 
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for details see Tables [7] and [9] 
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Table 10: Summary of the on-orbit Fermi LAT calibrations. 
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